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Abstract 

This note finds results for the intrabeam scattering growth rates 
for a bi-gaussian distribution. The bi-gaussian distribution is inter- 
esting for studying the possibility of using electron cooling in RHIC. 
Experiments and computer studies indicate that in the presence of 
electron cooling, the beam distribution changes so that it developes a 
strong core and a long tail which is not described well by a gaussian, 
but may be better described by a bi-gaussian. Being able to compute 
the effects of intrabeam scattering for a bi-gaussian distribution would 
be useful in computing the effects of electron cooling, which depend 
critically on the details of the intrabeam scattering. The calculation 
is done using the reformulation of intrabeam scattering theory given 
in [1] based on the treatments given by A. Piwinski [2] and J. Bjorken 
and S.K. Mtingwa [3]. The bi-gaussian distribution is defined below as 
the sum of two gaussians in the particle coordinates x,y, s,p x ,p y ,p s . 
The gaussian with the smaller dimensions produces most of the core 
of the beam, and the gaussian with the larger dimensions largely pro- 
duces the long tail of the beam. The final result for the growth rates 
are expressed as the sum of three terms which can be interperted re- 
spectively as the contribution to the growth rates due to the scattering 
of the particles in the first gaussian from themselves, the scattering 
of the particles in the second gaussian from themselves, and the scat- 
tering of the particles in the first gaussian from the particles in the 
second gaussian. 
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1 Introduction 

This note finds results for the intrabeam scattering growth rates for a bi- 
gaussian distribution. 

The bi-gaussian distribution is interesting for studying the possibility of 
using electron cooling in RHIC. Experiments and computer studies indicate 
that in the presence of electron cooling, the beam distribution changes so 
that it developes a strong core and a long tail which is not described well 
by a gaussian, but may be better described by a bi-gaussian. Being able 
to compute the effects of intrabeam scattering for a bi-gaussian distribution 
would be useful in computing the effects of electron cooling, which depend 
critically on the details of the intrabeam scattering. The calculation is done 
using the reformulation of intrabeam scattering theory given in [1] based on 
the treatments given by A. Piwinski [2] and by J. Bjorken and S. Mtingwa 
[3]. The bi-gaussian distribution is defined below as the sum of two gaussians 
in the particle coordinates x,y, s,p x ,p y ,p s . The gaussian with the smaller 
dimensions produces most of the core of the beam, and the gaussian with the 
larger dimensions largely produces the long tail of the beam. The final result 
for the growth rates are expressed as the sum of three terms which can be 
interperted respectively as the contribution to the growth rates due to the 
scattering of the particles in the first gaussian from themselves, the scattering 
of the particles in the second gaussian from themselves, and the scattering of 
the particles in the first gaussian from the particles in the second gaussian. 

2 Basic results for intrabeam scattering 

This section lists some general results which can be used to find growth rates 
for a beam with any particle distribution f(x.p). Following [3], growth rates 
will be computed for < p^pj > , where the <> indicates an average over all 
the particles in the bunch. From these one can compute the growth rates for 
the emittances, < 6j >. A result that holds in any coordinate system and for 
any particle distribution f(x.p) is given in [1] as 



5 < (piPj) > = N d 3 x 



Cij = TT 
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Wi = Pli+P2i (1) 

Nf(x,p) gives the number of particles in d 3 xd 3 p, where N is the number 
of particles in a bunch. 5 < {piPj) > is the change in < {piPj) > due to all 
particle collisions in the time interval dt. The invariants F(pi,p 2 ), A 2 , W 2 
are given by 

= c [(PiP 2 ) 2 -m 2 m 2 c 4 ] 1/2 
' ■m 1 m 2 c 2 

F( Pl ,P2) = 7i72c[(/?i - fy 2 - 0i x fy 2 ] 1/2 

A 2 = A 2 -A 2 , A = (E 1 -E 2 )/(2c) 

W 2 = W 2 - W 2 , W = + J B 2 )/c 

Eq.(l) is considerably simplified by going to the rest CS , which is the 
CS moving along with the bunch and the particle motion is non-relativistic, 
and putting a equal to the Coulomb cross section. One gets 

\ < 5(p u p lj ) > — N [ d 3 xd 3 p 1 d 3 p 2 f{x,p 1 )f(x,p 2 )2(3cC ij dt 
Po J 

Ai = -(Pli-P2i) 

(3c = \A\/m 

C i3 = 4(r /2^) 2 ln(l + (2^ 2 6 max /r ) 2 ) 
Po 

[|A|%-3A,A,] i,J = 1,3 
r = Z e jmc 

cot{e min /2) = 2p 2 b max /r (2) 

b m ax is the largest allowed impact parameter in the center of mass CS. It 
has been asumed that one can replace ln(l + (2(3 2 b max /r ) — 1 by ln(l + 
(2f3 2 b max /r ). 

In Eq.(l), the original 11-dimensional integral which arises from intra- 
beam scattering theory has been reduced in [1] to a 9-dimensional integral 
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by integrating over all possible scattering angles. In [1] this reduction was 
done for any particle distribution, f(x,p). In [3], Bjorken and Mtingwa first 
do the integration over x,p\,p2 using a simple gaussian distribution before 
doing the integration over the scattering angles and no general result for 
doing this reduction for any f(x,p) is given. In [2] Piwinski computes the 
growth rates for the emittances < > instead of for < piPj >. A general 
result for reducing the integral by integrating over all possible scattering 
angles, for any f(x,p), for the growth rates of < q > is given. However, 
using this result for a complicated distribution like the bi-gaussian would be 
difficult. 

3 Gaussian distribution 

We will first consider the case of a gaussian particle distribution. This will 
provide a more simple example of using the results in the reformulation given 
in [1] and of the methods used to evaluate the integrals. Afterwards, the same 
procedures will be applied to the case of the bi-gaussian distribution. 

Let Nf(x,p) gives the number of particles in d 3 xd 3 p, where N is the 
number of particles in a bunch. For a gaussian distribution, f(x,p) Is given 



by 



f(x,p) = -exp[-S(x,p)] 



r 




(3) 



S = 




S x — —£x(Xf3, Xp) 




Xp = x - D(p- p )/po 
x'p = x' - D'(p-p )/po x' = p x /p 
x') = [x 2 + {(3 x x' + a x x) 2 }/(3 x 



S, 



y 



—£y(y,y) y' = p y /po 



+ (p y y> + a y y) 2 }/% 



•y 
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S s 


= — e s (£ 




- s c , (p 


-Po)/Po) 






- s c , (p 


-Po)/Po) 


"S<- 




- s c , (p 


-Po)/po) 


= [(*-■ 






ft 


= o" s /o"p 








= 2a s a p 



(4) 

D is the horizontal dispersion. D' = dD/ds. A longitudinal emittance has 
been introduced so that the longitudinal motion and the transverse motions 
can be treated in a similar manner. s c locates the center of the bunch. 

T can now be computed using Eq.(l).This will provide an example how 
the integrals are done in this paper. The integration methods used here are 
somewhat more complicated than those used in [3] but they will also work 
for the more complicated bi-gaussian distribution. 

r = J d 3 xd 3 p exp[-S x - S y - S s ] 

Writing T as T = T y T xs and computing Y y first because this part is simpler, 





= J dydp y exp[-S y ] 


Sy 


= —t y {y,y') y' = p y /po 


e y (y,y') 


= [y 2 + 2 + (Pyy' + a y y) 2 }/Py 


Vy 


= y/ \[Kn Pw = (/W + a y y)l \[K 


dydpy 


= podrjydprjy 


r„ 


= Po j driydpryy exp[-{rf y + p 2 m ) / e y ] 




= neypo 



(5) 



Now for the remaining integral we have 

r xs = J dxdp x dsdp s exp[—S x — S s 



T xs = J dsdp s exp[-S s ] J dxdp x exp[-S x 
Make the transformation 

xp = x - D(p- p )/Po 

x'p = x' - D'(p-p )/p 

d = Px/PO, x'(3=P(3x/P0 

dxdp x = podxpdx'p 

Sx — —^x{.Xp,x'p) 



J dxdp x exp[— S x ] = TTe x po as in evaluating T y 
p ~ p s in the Lab. CS and 

T = ir 3 e s 6 x e y p 3 (6) 

4 Growth rates for a Gaussian distribution 

In the following,the growth rates are given in the Rest Coordinate System, 
which is the coordinate system moving along with the bunch. Growth rates 
are given for < p(Pj >. From these one can compute the growth rates for 
< €i >. Using the general result, Eq.(2), one gets 



\ < S(piPj) > = ^ f d 3 xd 3 p 1 d 3 p 2 exp[-S(x,p 1 ) - S(x,p 2 )]2f3cC ij dt 
Po * J 



1 

A = g (Pi "Pa) 

(3c = \A\/m 
lir 

C l3 = — (r /2(3 2 ) 2 \n(l + (2(3 2 b max /r ) 2 ) [\A\% - 3A.A,.] i,j 
Po 

7 2 2/ 2 

r = Z e jmc 
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r = 7r 3 e s e x e y p 3 (7) 
Transform to W, A 

pi = y + pz = y~ 

= Pi +p 2 , A = — - — 

d 3 Pl d 3 p 2 = d 3 Wd 3 A (8) 

We will first do the integral over d 3 x and over d 3 W. For the y part of the 
integral 

s y (y,piy) = —£ y (y,y'i), y'i=Ply/Po 

e y 

£ y (y,y'i) = [y 2 + (M + ^ y y) 2 }/p y 

Sy(y,Pi y ) = [y 2 + {(3 y {^ + A y )/p Q + a y y) 2 }/((3 y e y ) 

S y (y,p ly ) + S y (y,p 2y ) = (2y 2 /(3 y + 2((3 y (W y / Po )/2 + a y yf / f3 y 

+2(3 2 (A y /p ) 2 /(3 y ))/e y 

Make the transformation 

Vy = V2y/Jp~ y , p vy = V2({3 y (W y /p )/2 + a y y)/^P~ y 
dydW y = p dr] y dp m (9) 

Integrate over dy, dW y 

J dydW y exp[-S y (y,p ly ) - S y (y,p 2y )] = p J drj y dp m 

exp[-^ + ^ !!i; 



V'y + P\y + 2p 2 (A y / P0 ) 2 /P y ] 



■y 



2(3 

p ne y exp[ ^(A y /p ) 2 } 

e y 

p ire y sxp[-Ry\ 



Ry = 2A (\/Po) 2 (10) 
e y 
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In doing the remainder of the integral, the integral over dxdW x dsdW s we 
will do the integral over dxdW x first and then the integral over dsdW s . Note 
that the integral is being done in the Rest CS and in the expression for S x 
one has to replace p — po ~ p s — po in the Lab CS by 7p s in the Rest CS. 
Remember also that f(x,p) is an invariant (see [1]) One finds for S x (x,pi x ) 



S x (x,p lx ) = {{x- 1 DW s /2- 1 DA s ] 2 + [(3 x {W x /2 + A x - 1 D'W s /2- 1 D'A s ) + 
a x (x - >yDW a /2 - 1 DA s )] 2 }/(f3 x e x ) 

Wi = Wi/po A l = A/p 

S x (y,p lx ) = {[x- 1 DW s /2- 1 DA s } 2 + [p x (W x /2- 1 D'W s /2) + 
a x (x - -yDWs/2) + {(3 X A X - ^DA,)} 2 } / {^ x e x ) 
D = p x D' + a x D (11) 



we then find for 



S x (x,Pix) + S x {x,p 2x ) 



S x (x,p lx ) + S x (x,p 2x ) 



{2[x - -yDWs/2] 2 + 2 1 2 D 2 A 2 S + 

2[(3 x {W x /2 - jD'Ws/2) + a x {x - 1 DW s /2)] 2 + 

2[(3 x A x - 1 DA s ] 2 }/((3 x e x ) (12) 



Now make the transformations 



Vx 

dxdW x 



x 




(13) 



Doing the integral over dxdW x one finds 
/ dxdW x exp[-S x (x,p lx ) - S x {x,p 2x )} -■ 



Po J dr] x dp VxX 
exp[-{vl + P 2 Vx + 

2[ 1 2 D 2 A 2 s + (p x A x 

p ire x exp[-R x ] 



7J DA s ) 2 ]/^)}/ e -] 



Rx — 



2[ 1 2 D 2 A 2 S + (p x A : 



7DA s ) 2 ]/(p x e x ) 



(14) 
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Now do the integral over dsdW s . One may note that the form of the 
intgral here is similar to the integral done over dydW y . The result is then 
the same with the proper sustitutions of s for y. 

J dxdW s exp[-S s (s,pi s ) - S s (s,p 2s )] = p ire s exp[-R s ] 

R s = ^£(A s / Po ) 2 (15) 

Note that the term f3 s ((p — Po)/po) 2 in S s in the Lab. CS has to be replaced 
by l 2 Ps(Ps/Po) 2 in the Rest CS. 

Using Eq.(7), one gets the result for the growth rates in the Rest CS for 
a gauusian distribution. 



1 d JV / r 

-$— + < PiPj > = y J exp[-R\Cij 

2n 

Cij = — (r /2/3 2 ) 2 (|A|% - 3^)2/^ ln[l + (2^ 2 b max /r f 
Po 



pi dt 



P = A)7o|A/p, 



o 



r = Z 2 e 2 /Mc 2 
T = 7r 3 e~e~e~Po 

R = R x + R y + R s 

R x = ^h 2 D 2 A 2 s + (p x A x -^DA s ) 2 ]/pl 

Px^x 



D = (3 x D' + a x D 

Ry = ^K/Pl 
e y 

R s = 2 -^i 2 A 2 Jp 2 (16) 



The integral over d 3 A is an integral over all possible values of the relative 
momemtum for any two particles in a bunch. [3q, 70 are the beta and gamma 
corresponding to po, the central momemtum of the bunch in the Laboratory 
Coordinate System. 7 = 70 
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The above 3-dimensional integral can be reduced to a 2-dimensional inte- 
gral by integrating over |A| and using d 3 A = \A\ 2 d\A\sin9d9d(p. This gives 



p\dt 



2 ^ < PiPj > 



93 

9i 

92 

C 



r 



N , 



c 

cosO = g s 
sin9cos(f> = g x 
sin9sin(p = g y 
^lf3lb max /r 



2/5 7oc / sinddOdcj) - Sgigj) 



F = 
F = 

F„ = 



-^-\l 2 D 2 gl^^ x g x - 1 Dg a f\ 



-Psl 2 9 



2.2 



In obtaining the above, one uses z = \ A\ 2 ,dz = 2|A|d|A| and 



1 



dz exp[—Fz]ln[Cz] = —[In 
o F 



- .5772] 

For Z = 80,A = 200, 7 = 100, b max = lcm, log 10 C = 18.6 



(17) 



5 Bi-Gaussian distribution 

The bi-gaussian distribution will be assumed to have the form given by the 
following. 

Nf(x,p) gives the number of particles in d 3 xd 3 p, where N is the number 
of particles in a bunch. For a bi-gaussian distribution, f(x,p) Is given by 



f(x,p) = ^^exp[-S a (x,p)] + ^^-exp[-5 6 (x,p)] 



10 



1 a 7T ^sa^xa^yaPo 

= 7r 3 e~ 6 e~ 6 e~ 6 po ( 18 ) 

In the first gaussian,to find T a , 5 a then in the expressions for T, S, given 
above for the gaussian distribution, replace e x , e y , e s by e xa , e^ , e~ sa . In the 
second gaussian, in the expressions for T, S, replace e x , e y , e s by e x b, e^, e^,. 
In addition. N a + N b = N. This bi-gaussian has 7 parameters instead of the 
three parameters of a gaussian. 



6 Growth rates for a Bi- Gaussian distribu- 
tion 

In the following,the growth rates are given in the Rest Coordinate System, 
which is the coordinate system moving along with the bunch. Growth rates 
are given for < piPj >. From these one can compute the growth rates for 
< €i >. Starting with Eq.2 and using the f(x,p) from Eq.18, one gets 



4 < 6 (PiPj) > 

Po 



Nn 1 



N h 1 



A 

/3c 



a 



r 



J d 3 xd 3 p 1 d 3 p 2 ^^exp[-S a (x,pi)} + ^^expl-Sbix^p!)} 

2(3c Cij dt 

\{p"i -P2) 
\A\/m 

^(r /2p 2 ) 2 \n(l + (2p 2 b max /r ) 2 ) [\A\% - 3A,A,-] i,j = 1,3 
Po 

Z 2 e 2 /mc 2 (19) 



The term in the integrand which contains exp[— S a (x , pi) — S a (x , P2)] is similar 
to the integrand for the gaussian distribution except that q are replaced by 
e ia and leads to the same result as that given by Eq.(16) for the gaussian 
beam except that R has to be replaced by R a where R a is obtained from R 
by replacing q by e ia . The term containing exp[—Sb(x,pi) — Sb(x,p2)] can be 
evaluated in the same way leading to the same result as that given by Eq.(16) 
for the gaussian beam except that R has to be replaced by Rb where Rb is 
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obtained from R by replacing e$ by e^. The only terms that need further 
evaluation are the the two cross product terms. The two cross product terms 
are equal because of the symmetry of p\ and p 2 in the rest of the integrand. 
This leads to the remaining integral to be evaluated 

J d 3 xd 3 p 1 d 3 p 2 2N ^ b YTjrexpl-Saix, Pi) - S b (x,p 2 )} 2(3c dj 

In evaluating this integral, we will use the same procedure as was used 
for the gaussian distribution. We will first transform to W, A from pi,p 2 (see 
Eq.(8). We will then do the integral over d 3 x and over d 3 W. For the y part 
of the integral one finds , 

S ya (y,Piy) = {y 2 + [(3 y (W y /2 + A y ) + a y y] 2 } / ((3 y e ya ) 

Wy = W y /P Ay = Ay/ P0 (20) 

One then finds that 



S ya (y,Piy) + Sy b (y,p 2 y) = {2y 2 / (3 y + 2[(3y{W y /2 + a y y)\ 2 /(3 y 

+2p 2 y Al/p y )}/e yc + 
{A(P y A y )(P y W y /2 + a y y)/p y }/e yd 



1 


1 


, 1 


+f> 


e yc 


= 2 1 


y a 


1 


1 


, i 


-±) 


tyd 


= 2 1 


e ya 


tyb 



(21) 



Make the transformation 

Vy = V2y/\/Ty, 
dydWy = podi] y dp vy 



Pvy 

V2(f3 y Wy/2 + a y y)/M 



(22) 



Integrate over dy, dW y 

J dydW y ecp[-S ya (y,p ly ) - S yb (y,p 2y )] = p J drj y dp ri 



exp[- 



V 2 y+P 2 m + W 2 y (\/P0) 2 /Py 



+ 



-yc 
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p 7ie yc exp[-—A 2 + - 



= Poire yc exp[-R yc + i?^] 

p _ 2 A/a2 
itj/c - — A y 

Ry d = J^t^l (23) 

tydl t yc 

The exponent R yc — i?^ has to be positive. This can be made more obvious 
by noting that 

1 1 1 

-2 —2 — — 

e yc e yd e ya e yb 

In doing the remainder of the integral, the integral over dxdW x dsdW s we 
will do the integral over dxdW x first and then the integral over dsdW s . Note 
that the integral is being done in the Rest CS and in the expression for S x 
one has to replace p — p ~ p s — p Q in the Lab. CS by 7p s in the Rest CS. 
Remember also that f(x,p) is an invariant (see [1]) One finds for S xa (x,pi x ) 

S xa (x,p lx ) = {[x- 1 DW s /2- 1 DA s ] 2 + [p x (W x /2 + A x - 1 D'W s /2- 1 D'A s ) + 
a x (x - 1 DW s /2 - 1 DA s )] 2 }/(p x e x a) 

Wi = W./po A t = A/ Po 

S xa (y,Pi x ) = {[x- 1 DW s /2- 1 DA s ] 2 + [(3 x (W x /2- 1 D'W s /2) + 
a x (x - -yDWs/2) + (f3 x A x - 1 DA s )] 2 }/((3 x e xa ) 
D = p x D' + a x D (24) 

we then find for S xa (x,p lx ) + S xb (x,p 2x ) 

S xa (x,p lx ) + S xb (x,p 2x ) = {2[x- 1 DW s /2] 2 + 2 1 2 D 2 A 2 s + 

2[p x (W x /2 - ^D'Ws/2) + a x (x - ^DWs/2)} 2 + 



2[p x A x - 1 DA s } 2 }/((3 x e xc ) + 
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{-4 7 £A s [x- 7 W s /2]/& 

+A{(3 X A X - 1 DA s )\(3 x {W x /2 - ^D'W s /2) + 

a x (x - 1 DW s /2)]/p x }/e xd 



x* = V2x -^DW S /V2 p* x = W x /V2- 1 D'W s /V2 



Vx = x*/Jp x p VxX = (p x p* x + a x x*)/Jp x 



Doing the integral over dxdW x one finds 



dxdW x exp[-S xa (x,p lx ) - S xb (x,p 2x )] = p dr} x dp VxX 



exp[-{rjl + p 2 Vx 

+2[ 1 2 D 2 A 2 + (p x A x - 1 DA S ) 2 ]/{3 X }/^ 

+{-4 7 DA s ^//% 

+4{(3 X A X - 1 DA s )p VxX /^f x }/e xd } 

= p ne xc exp[-R xc + R xd ] 

R xc = 2[ 1 2 D 2 A 2 s + ^ x A x - 1 DA s ) 2 ]/((3 x e xc ) 
R xd = 2{[- 1 DA s ] 2 

+[(p x A x - 1 DA s )] 2 } 



Now do the integral over dsdW s . One may note that the form of the 
intgral here is similar to the integral done over dydW y . The result is then 
the same with the proper sustitutions of s for y. 

/ dxdW s exp[-S sa (s,p ls ) - S sb (s,p 2s )} = p ire sc exp[-R sc + R sd ] 



(25) 



Now make the transformations 



dxdW, 




(26) 



(27) 
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Note that the term (3 s ((p — po)/po) 2 in S s in the Lab. CS has to be replaced 
by l 2 Ps(Ps/po) 2 in the Rest CS. 

Putting all the above results, for the bi-gaussian distribution, together 
one gets the final result 



2n 

Cij = — (r /2^ 2 ) 2 (|A|%--3A,A,)2^c ln[l + (2p 2 b max /r ) 2 } 
Po 

P = /5o7o|A/p | 
r = Z 2 e 2 /Mc 2 

1 1/1 1 \ . 

— = 9 ( — + — J i = x,y,s 

JL = 

r = Z 2 e 2 /Mc 2 

Ta 7T tsa^xa^yaPo 



SCI 



^ = -r^— [7 2 -D 2 A 2 + (/3 X A X — ^fDA s ) 2 }/pl 

Px^-xa 

D = p x D' + a x D 

I y y a 

Rsa = -^-&VA 2 / P 2 

Ps^sa 
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T = 

T r . = 



Rxc Rxd 



Ty R 

T s = R x 



yc 



Ryd 
Rsd 



R xd = 2{[- 1 DA s ] 2 

+ [(P X A X - 1 DA S )] 2 } 

I (Px^xd/^xc) 

_ 2[3y ~ 2 

tydl t yc 
R.sd = _n ^7Z A, 



e 



1d/ e s. 



A, = A,/p (29) 

i? a , _R C are each the same as _R a given above except that e7 a are replaced 
by €i a , elb, e ic respectively. The same remarks apply to r a , F b , T c 

The above 3-dimensional integral can be reduced to a 2-dimensional inte- 
gral by integrating over |A| and using <i 3 A = \A\ 2 d\A\sin6d6d(f). This gives 



Id ( tq \ f 

— -f. < PiPj > = 2 ^Po Trhv 2 PoloC / sin9d9d(j) (5„ - 3^-) 
p at YZ7 Po/ J 

2 



N a N b r c 1 c 

93 = cosO = g s 

gi = sinOcoscj) = g x 

g 2 = sin9siri(p = g y 

C = 2^lf3lb max /r Q 

Fi = R t /(\A\/ Po ) 2 i = a,b,c 
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G = T/(\A\/ Po ) 2 



(30) 



F a , Ff,, F c are each the same F that was defined for the Gaussian distri- 
bution except that the q are replaced by C ia , e7&, e ic respectively. 

The above results for the growth rates for a bi-gaussian distribution are 
expressed as an integral which contains 3 terms, each of which is similar 
to the one term in the results for the gaussian distribution. These three 
terms may be given a simple interpertation. The first term represents the 
contribution to the growth rates due to the scattering of the N a particles of 
the first gaussian from themselves, the seond term the contribution due to 
the scattering of the iVj, particles of the second gaussian from themselves, 
and the third term the contribution due to the scattering of the N a particles 
of the first gaussian from the Nb partcles of the second gaussian. 



7 Emittance growth rates 

One can compute growth rates for the average emittances, < e« > in the Lab- 
oratory Coordinate System, from the growth rates for < piPj > in the Rest 
Coordinate System. In the following , dt is the time interval in the Laboratory 
System and dt is the time interval in the Rest System, dt = jdi 



d B x d 9 . 9 D 2 + D 2 d 9; n n ~ d . 9 

j t <£s> = Ps-yi < vllvl > (31) 

To derive the above results, the simplest case to treat is that of the vertical 
emittance. The verical emmitance is given by 

e y (y,y') = [y 2 + (f3 y y' + a y y) 2 }/ (3 y 
5e y = (3 y 5(y' 2 ) 



(32) 
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In Eq.(32), y' = p y /po, Se y is the change in e y in a scattering event. 
For the longitudinal emittance one finds 

e s = [s 2 /l 2 + (PslPs/Po) 2 }/Ps 
5e s = Ps^iiPs/Po) 2 

j t <£s> = Pslj- < pi/pi > 



(33) 



In Eq.(33), s,p s are the coordinates in the rest system and I have used the 
relationship {p-po) L AB = {iPs)rest 
For the horizontal emittance one finds 

e x = {[x-'jDps/po] 2 + [p x (p x /p -'jD'p s /po) + a x (x-'jDp s /p )] 2 }/p x 
= {[x- 'jDps/po] 2 + [P x p x /p + a x x - DiPs/pq] 2 }/ 'fix 
= {x 2 + {jDp s /pq) 2 - 2x>yDp s /p + ((3 x p x /p + a x x) 2 + (D>yp s /p ) 2 - 
2(p x p x /po + a x x)(D~fp s /p )}/p x 
5e x = 5{p 2 x (p x /p ) 2 + 1 2 (D 2 + D 2 )(p s /p ) 2 -2P x D 1 (p x /p )(p s /p )}/p x 
d B x d 9;9 D 2 + D 2 d 9;9 ~d . 9 

It < € ' > = -Jt < pJp " > +— — 7 S < R/Po > ~ 2D 7i < p ' p ' /p ° > 

(34) 

In the result for 5e x , the terms that are linear in p x or p s have been dropped 
as they do not contribute to < 5e x > . In a scattering event involving two 
particles , the 5p x of one particle is equal and opposite to the 5p x of the other 
particle. This is also true for p s . 
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